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A POWERFUL HYDRODYNAMIC BOOSTER FOR RELATIVISTIC JETS 
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ABSTRACT 

Velocities close to the speed of light are a robust observational property of the jets observed in microquasars 
and AGNs, and are expected to be behind much of the phenomenology of GRBs. Yet, the mechanism boosting 
relativistic jets to such large Lorentz factors is still essentially unknown. Building on recent general-relativistic, 
multidimensional simulations of progenitors of short GRBs, we discuss a new effect in relativistic hydrodynamics 
which can act as an efficient booster in jets. This effect is purely hydrodynamical and occurs when large velocities 
tangential to a discontinuity are present in the flow, yielding Lorentz factors V ~ 10 2 — 10 3 or larger in flows with 
moderate initial Lorentz factors. Although without a Newtonian counterpart, this effect can be explained easily 
through the most elementary hydrodynamical flow: i.e. a relativistic Riemann problem. 

Subject headings: galaxies: jets — gamma rays: bursts — hydrodynamics — relativity — shock waves 



pointed out by Rezzolla & Zanotti (2002), and understanding 
its influence on the complex dynamics of a relativistic jet can 
be very difficult, we resort to the simplest hydrodynamical flow 
in which this effect can occur: i.e., a Riemann problem. 

2. THE RELATIVISTIC BOOSTING MECHANISM 

We recall that a Riemann problem consists of determining 
the evolution of a fluid which is composed at some initial time 
by two uniform states (i.e., a "left" and a "right" state) with dif- 
ferent and discontinuous hydrodynamical properties: the rest- 
mass density p, the pressure p, the specific internal energy e, 
the specific enthalpy h = 1 + e/c 2 + pj (pc 2 ) and the compo- 
nents of the velocity normal v n and tangential v l to the initial 
discontinuity. Using this initial setup and the equations of rel- 
ativistic hydrodynamics, Rezzolla et al. (2003) revealed that, 
in contrast with what happens in Newtonian hydrodynamics, a 
smooth transition from one wave-pattern to another can be pro- 
duced by simply varying the velocities which are tangential to 
the initial discontinuity u*, while maintaining the initial states 
unmodified. More specifically, for each set of initial states there 
exists a critical value for the initial tangential velocities at which 
a shock is transformed into a rarefaction wave or viceversa. 

The discovery of this relativistic effect is not purely academic 
and without direct physical implications. Indeed, an impor- 
tant consequence of this change of waves is the possibility of 
producing a large amplification of the local the Lorentz factor 
T = [1 — (v 1 ) 2 — (V 1 ) 2 ] -1 / 2 when, past a critical value for the 
tangential velocity, a shock wave transforms into a rarefaction 
wave. To illustrate how this can happen in practice it useful 
to consider a Riemann problem with initial states that would 
mimic the conditions for the generation of a GRB or the prop- 
agation of a plasma-blob along a relativistic jet (either extra- 
galactic or in a microquasar). This can be done in practice by 
considering a left state which is much hotter, with larger pres- 
sure and smaller rest-mass density than the right state. Stated 
differently, this corresponds to a left state having a much larger 
specific enthalpy than the right one. Such a choice is justi- 
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1. INTRODUCTION 

Determining the mechanism accelerating high-energy 
plasma to the relativistic velocities observed in jets from Ac- 
tive Galactic Nuclei (AGNs) and microquasars, and expected in 
gamma-ray bursts (GRBs) remains a problem whose solution 
has long been sought in any consistent modelling of these astro- 
physical objects. We here show how recent general-relativistic 
simulations and novel effects of relativistic hydrodynamics can 
provide important clues in this search. 

In crude but not oversimplified terms, the dynamics of a rel- 
ativistic hydrodynamic jet in an external medium can be assim- 
ilated to the motion of two fluids, one of which is much hot- 
ter and at higher (or equal) pressure than the other one, and is 
moving with a large tangential velocity with respect to the cold, 
slowly moving fluid. In this scenario, a blob of hot plasma 
moving outwards along the jet generically produces a hydrody- 
namical structure, in the direction perpendicular to the motion, 
composed of a "forward" shock (moving away from the jet- 
axis) and of a "reverse" shock (moving towards the jet-axis). 
This pattern of waves, which we indicate as ^_5C6U, where 
refers to the reverse shock, 6U to the forward one and C 
to the contact discontinuity between the two, is routinely pro- 
duced in numerical simulations of relativistic jets. Surprisingly 
however, because unexpected from Newtonian hydrodynamics, 
this pattern can change for sufficiently large tangential veloci- 
ties and in this case the inward-moving shock is replaced by a 
rarefaction wave thus producing a pattern. Convincing 

evidence that this process is likely to occur under realistic con- 
ditions is now emerging from special relativistic calculations 
of extragalactic jets (Aloy et al. 2003), as well as from gen- 
eral relativistic, multidimensional simulations of progenitors of 
short GRBs (Aloy et al. 2005). The latter case is schematically 
shown in Fig. 1 which shows the flow-structure produced in nu- 
merical simulations involving a thick accretion torus around a 
black hole and in which a ultrarelativistic outflow is generated. 

Because disentangling this change in the wave-pattern, first 
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fied by the fact that the jet is expected to be very underdense 
with respect to the external medium (which could be given, for 
instance, by a dense accretion torus) and to be overpressured 
because a substantial change in the jet radius must happen be- 
tween its launching site and the location at which it is detected 
z. For GRBs, for instance, one expects z w 10 13 — 10 14 cm, 
while z Kt 10 16 — 10 17 cm for microquasars and z w 10 18 cm 
for extragalactic jets. 




FIG. 1 . — Schematic illustration of the flow structure in the acceleration of a 
short GRB produced by an accretion torus orbiting around a stellar-mass black 
hole. The arrows mark the direction of fluid velocity at the rarefaction head 
(yellow paraboloidal surface), indicating that collimation of the fluid tends to- 
wards the black-hole rotation axis (black arrow). A large boost is produced 
in the region between the rarefaction head and the contact discontinuity (red 
surface) separating the relativistic outflow from the shocked external medium. 

Keeping these initial states fixed, we can then use the tan- 
gential velocity in the left state u* as a free parameter, whose 
tuning allows to go from one wave-pattern to the other. More 
specifically, let us consider the right state to be a cold fluid with 
a large rest-mass density and essentially at rest, i.e., with initial 
hydrodynamical properties: p n w 10~ 6 /9 ext c 2 , p R « W~ 2 p mt , 
wO, k 0, where p oxt is a normalization constant which 
allows us to be scale-free. For convenience, herafter we will 
assume c = 1 and provide explicit values for p cxt when dis- 
cussing the astrophysical implications in Sect. 3. The left state 
will have comparable but larger pressure, will be underdense, 
with a small initial tangential velocity and a normal one close 
to the speed of light, i.e., p L > p Tt , p L <C p R , u" 0. 

The solution of the Riemann problem produced by these two 
states after an arbitrary time is shown in the left panel of Fig. 2 
(We recall that the solution of the Riemann problem is self- 
similar, i.e., the structure of the solution in the different waves 
does not depend on time.). The solution is of type ^SC6U and 
thus consisting of a shock propagating to the left, of a shock 
propagating towards the right and of a contact discontinuity 
at x w 0.369 separating the two waves. This double-shock 
structure has a well-known Newtonian counterpart and is com- 
monly observed in numerical simulations when there is an al- 
most pressure-matching between the jet and the environment. 
Note that the acceleration at the transition layer between the jet 



and the external medium is very small and that the Lorentz fac- 
tor reached in the state between the left-propagating shock and 
the contact discontinuity, T* , is smaller than the initial Lorentz 
factor in the left state T L . This is indicated as (T* ) max in Fig. 2 
and reaches a value (T*) max 2.943 against an initial one in 
the left state T L = 20. 

While keeping the right state to be the same, let us now con- 
sider the jet to be considerably over-pressured and with almost 
no velocity component perpendicular to the jet direction. The 
left state can then be parameterized as: p L p R , p L Cp B , 
u"«0, «1, and the solution of the Riemann problem in 
this case is shown in the right panel of Fig. 2. This is now of 
type ^JZCS-t and consists of a rarefaction wave propagating to 
the left, of a contact discontinuity at x w 0.027 and of a shock 
propagating towards the right. This change in wave-pattern is 
purely relativistic and it is accompanied by a large increase in 
the Lorentz factor across the rarefaction wave, with the maxi- 
mum value (r* ) max ps 10 3 reached at the contact discontinuity. 

This is a remarkable feature of relativistic hydrodynamics: 
namely, the production of a very strong rarefaction near the con- 
tact discontinuity and the combined conservation across such a 
wave of the specific enthalpy and of the Lorentz factor (i.e., of 
the quantity hT), has the consequence of accelerating the fluid 
to ultrarelativistic regimes with T ~ 10 4 or larger. In essence, 
this boosting mechanism efficiently converts into kinetic energy 
the work done on the jet by the external medium. 

It is worth stressing again that the occurrence of these large 
fluid accelerations is not the result of peculiar initial conditions 
in a simplified hydrodynamical flow. Rather, accelerations pro- 
duced through this effect are a robust feature in numerical sim- 
ulations of relativistic jets (Aloy et al. 2003, 2005). In par- 
ticular, the recent simulations of Aloy et al. (2005) show that, 
in spite of the complex multidimensional dynamics produced 
by the generation and propagation of the jet, the basic feature 
of this relativistic effect is preserved: i.e., the growth of the 
Lorentz factor close to the contact discontinuity separating a 
high-enthalpy and high-velocity state from a low-enthalpy one. 

As an illustrative example, we report in Fig. 3 results from 
an axisymmetric simulation in which the ultrarelativistic jet is 
overpressured with respect to the external medium and its ve- 
locity is almost parallel to the jet/external medium interface 
(Aloy et al. 2005). Shown along the polar ^-direction at a dis- 
tance of 3.2 x 10 7 cm from the black hole are the logarithm 
of the rest-mass density (dashed line) and the profile of the 
Lorentz factor (solid line). The latter, in particular, reaches a 
maximum T w 29 in the left-going rarefaction wave moving 
towards the axis. These simulations have also revealed that al- 
though the interaction between the jet and the external medium 
is prone to Kelvin-Helmoltz and shear driven instabilities (Aloy 
et al. 2002), the growth of the Lorentz factor imprinted during 
the initial phases of the jet propagation, i.e., when the jet in- 
teracts with the lateral borders of the accretion torus, remains 
preserved in the subsequent evolution of the flow. 

Let us now illustrate in detail how the final fluid accelera- 
tion depends on the initial conditions by assuming that the right 
state is held fixed. In this way, in fact, the transition from one 
wave-pattern to the other and the boost in the flow depend on 
two distinct but related factors: the value of the pressure in the 
left state p L and the value of the tangential velocity there, t>* . 

In particular, it is possible to show that for very small left- 
state pressures (i.e., p L < 10~ 6 for the initial states used 
above) the system evolves to produce a Riemann structure of 
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FIG. 2. — Left panel: Profiles for the Lorentz factor (left) and rest-mass density (right) resulting from a prototypical Riemann problem yielding a ^SCS-^ 
pattern. The initial states are (p L , p L , v™ , T L ) = (10~ 3 , 10~ 4 , 0.99, 20) and (p R , p R , v™ , T R ) = (10~ 6 , 10~ 2 , 0, 1). Note that no boost is produced 
at the contact discontinuity (i.e., (r* ) max < T L ). Right panel: Same as in the left panel but for a prototype ^JZCS-^ wave-pattern. The initial states are 
(p L ,p L ,u£,r L ) = (10 3 , 10- 4 ,0,20) and (p R , p R , v™, T R ) = (10~ 6 , 10~ 2 , 0, 1). Note in the left panel the very large boost (i.e., (r*) max > r L )which 
peaks at the contact discontinuity at x fs 0.027. 



the type ^SCSL, also when the normal velocity is zero {i.e., 
v™ w 0). On the other hand, for very large left-state pres- 
sures (i.e., p L > 1 for initial states used above), the system 
evolves to produce a Riemann structure of the type ^_7£C<S_> 
also when the normal velocity is ultrarelativistic (i.e., v" w 1). 
For all the values of p L within this range, a critical normal ve- 
locity (f™) c exists below which the left-propagating shock is 
replaced by a rarefaction wave and a large boost takes place. 
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FIG. 3. — Example of the growth of the Lorentz from a multidimensional 
simulations of ultrarelativistic jets generated in post-neutron star mergers (data 
from model A09 of Aloy et al. 2005). This is to be compared with the similar 
results shown in the right panel of Fig. 2. 

All of this is summarized in Fig. 4, where the maximum val- 
ues of the Lorentz factors are shown as function of the initial 



value of the normal velocity in the left state u" (Note that the 
value of T L is held fixed so that the increase of u" corresponds 
to a decrease of u* .). Different curves refer to different values 
of the pressure in the left state p L , with solid lines indicating 
solutions of the Riemann problem of the type ^JZCS-^, dashed 
lines to solutions of the type <__SCS->, and filled circles mark- 
ing the critical normal velocity (t>") c at which one wave-pattern 
transforms into the other. 

3. ASTROPHYSICAL IMPLICATIONS 

A first important astrophysical implication of the relativistic 
effects discussed here is that of requiring less extreme energy 
contents at the launching sites of relativistic jets and GRBs. 
This hydrodynamic mechanism, in fact, can boost a fluid to ul- 
trarelativistic regimes very efficiently and thus the Lorentz fac- 
tors either expected or observed in astrophysical sources can 
be explained with initial "left-states" that have considerably 
smaller pressures and are thus easier to produce in practice (in 
case the jet acceleration is mostly driven by pressure gradients). 

To be more specific, for an ultrarelativistic jet expected in 
a short GRB, the conditions close to the generation site (i.e., 
above the poles of a rapidly rotating black hole) could be as- 
sociated to a left state with p L ~ (10~ 2 — 10)p ext , i>™ « 0, 
u* w 1 and p L /p R <C 1 [For a short GRB progenitor p cxt ~ 
1 - 10 3 g/cm 3 (Janka et al. 1999)]. In this case, a Lorentz 
factor as small as T L ~ 10 is sufficient to produce, through the 
acceleration across the rarefaction wave, the inferred Lorentz 
factors of r o b s > 100. Furthermore, because the pressure in 
the jet decreases more rapidly than the pressure in the external 
medium (i.e., p L oc z~ 4 , p R oc z~ 3 , where z is the distance 
from the black hole) the boosting effect will be confined to a 
region within which p L /p R >L For a black hole with a mass 
of ~ 3Mq this may happen for z < 300 km. We also note 
that if the evolution of these boundary regions is maintained 
until the afterglow phase, then some differences in the emission 
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properties are to be expected with respect to the standard model 
of the "universal" jet-profile (i.e., a top-hat profile with sharp 
edges) or of the "structured-jet" model (i.e., where a progres- 
sive decay of the energy and of the Lorentz factor is expected 
from the jet core to the external medium). The impact of these 
structural changes within the jet on the light curve, whose cal- 
culation goes beyond the scope of this letter, will help to collect 
evidence on the occurrence of this relativistic acceleration. 
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FIG. 4. — Maximum Lorentz produced in the Riemann problem (r* ) max as 
function of the normal velocity in the left state v™. The right state is held 
fixed and is given by (p R , p R , v™ , T R ) = (ltT 6 , 1CT 2 , 0, 1). The left 
state has fixed values in the rest-mass density and Lorentz factor (p L ,T L ) = 
(1(T 4 ,20), while the pressure p L is varied for each curve as indicated by the 
different labels. Continuous lines refer to a <_7J.CS-> pattern, while dashed 
lines to a <— SC<S— filled circles mark the critical normal velocity (*>™) c . 

In extragalactic relativistic jets, on the other hand, the 
Lorentz factors are more moderate, with r o b s < 20 (Vermeulen 
& Cohen 1994) and the physical conditions could then be as- 
sociated to a left state withp L ~ (1(T 6 - 10~ 4 )p ext , u™ « 0, 
u* « 1 and Pl/Pj, <~ 10~ 2 - 10~ 5 [For an extragalactic jet 
Poxt ~ 10~ 27 - 10~ 24 g/cm 3 (Ferrari 1998).]. In this case, a 
Lorentz factor T L ~ 1.5 is sufficient to produce the observed 
boosts. In addition, if a continuous channel of plasma is present 
together with to the detected superluminal ejections, the phys- 
ical conditions in galactic microquasar jets might be similar to 
the ones found in extragalactic jets, but with an even smaller 
inferred value of the Lorentz factor (r o b s < 5; Mirabel & 
Rodriguez 1999). In this case, values as small as T L ~ 1.2 
would be sufficient to yield the observed apparent speeds. 

A second important astrophysical implication of our results 
is of pointing out that, under suitable but realistic conditions, 



ultrarelativistic regimes can be reached very close to the jet 
launching site. Furthermore, while the boost originates and is 
larger in a small region at the interface between the jet and the 
external material, it is not restricted to a thin boundary layer 
but can rapidly involve the whole jet. This is because the left- 
going rarefaction wave can sweep up the jet material at rela- 
tivistic speeds and thus rapidly accelerate a considerable mass 
fraction in the jet. As an example, on the basis of the values de- 
duced from the Riemann problem considered in the right panel 
of Fig. 2, more than 10% of the total (baryonic) mass is acceler- 
ated to T > T L in less that 1/5 of the lateral jet-crossing time. In 
contrast, in other thermal acceleration mechanisms that also ex- 
ploit the existence of h ^> 1 (e.g., the fireball model for GRBs) 
some time is needed to speed-up the flow, with ultrarelativistic 
regimes being reached far from the launching site. 

It is important to underline that the boosting mechanism pro- 
posed here should not be considered an alternative to the fireball 
model. Rather, it complements it by providing a more detailed 
description of the relativistic flows at the launching sites. In- 
deed, we expect both processes of hydrodynamic acceleration 
to operate simultaneously, with a longitudinal acceleration be- 
ing produced by the expansion of a hot gas, and a the lateral 
one being produced by the relativistic effects discussed here. 

Interestingly, this boosting mechanism can also enhance par- 
ticle acceleration when a turbulent velocity shear layer devel- 
ops along a parsec-scale jet. These layers, in fact, whose ex- 
istence has been pointed out both observationally (e.g., Owen 
et al. 1989; Swain et al. 1998) and through numerical simula- 
tions (Aloy et al. 1999, 2000), represent natural sites for par- 
ticle acceleration, providing high-energy cosmic rays and in- 
fluencing the dynamics of relativistic jets in extragalactic radio 
sources by forming cosmic -ray cocoons (Ostrowski 2000). The 
efficiency of the acceleration process in these turbulent shear 
layers depends on the particle mean-free-path and on the veloc- 
ity structure (Stawarz & Ostrowski 2002). Our results suggest 
that particle production can be further amplified if the simple 
flow structure usually assumed in these calculations and similar 
to the one produced in a <_££!£_> is instead replaced by a more 
realistic one, such as the one occurring in a <JJZCS-^. 

As a final remark, we point out that the large hydrodynamic 
boosts reported here may be detected with laboratory exper- 
iments involving heavy-ion collisions. In these experiments, 
heavy ions are accelerated to ultrarelativistic velocities and col- 
lided. If the two beams are chosen to have different specific 
enthalpies and to collide with a nonzero impact parameter, an 
acceleration of the interacting layer could be produced as a re- 
sult of the relativistic boosting discussed here. 

MAA is a Ramon y Cajal Fellow of the Spanish MEC. Sup- 
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swellenastronomie" of the DFG. 
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